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Abstract. The fine calibration of the ISO-SWS detectors (Infrared Space Observatory - Short Wavelength 
Spectrometer) has proven to be a delicate problem. We therefore present a detailed spectroscopic study in the 2.38 
- 12 fim wavelength range of a sample of 16 AO - M2 stars used for the calibration of ISO-SWS. By investigating 
the discrepancies between the ISO-SWS data of these sources, the theoretical predictions of their spectra, the 
high-resolution FTS-KP (Kitt Peak) spectrum of a Boo and the solar FTS-ATMOS (Atmospheric Trace Molecule 
Spectroscopy) spectrum, both calibration problems and problems in computing the theoretical models and the 
synthetic spectra are revealed. The underlying reasons for these problems are sought for and the impact on the 
further calibration of ISO-SWS and on the theoretical modelling is discussed extensively. 
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Appendix A: Comparison between different 
observational and synthetic data (coloured 
plots) 



In this section, Fig. ?? - Fig. ??, Fig. ?? - Fig. ?? of 
the accompanying article are plotted in colour in order to 
better distinguish the different observational or synthetic 
spectra. 
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Fig. A.l. Comparison between the ISO-SWS data of a 
Cen A (black) and the synthetic spectrum (red) with steL 
lar parameters Teff = 5830 K, log g = 4.35, M = 1.3Mq, 
[Fe/H] = 0.25, ^ = LOkms-\ ^"^C/^^C = 89, £(C) = 
8.74, e(N) = 8.26, e(0) = 9.13 and Od = 8.80 mas. Some 
of the most prominent discrepancies between these two 
spectra are indicated by an arrow. 
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Fig. A. 2. Comparison between band 1 and band 2 of 
the ISO-SWS data of 7 Dra (black) and the synthetic 
spectrum (red) with stellar parameters Teft = 3960 K, log 
g = 1.30, M = 1.7M0, [Fe/H] = 0.00, ^ = 2.0km s-\ 
12C/13C = 10, £(C) 8.15, £(N) = 8.26, e(0) = 8.93 and 
0d = 9.98 mas. 
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Fig. A. 4. Comparison between the ISO-SWS observation 
and the synthetic spectrum of 7 Dra at a resolving power 
of 1500 (bein g the most conser vative theoretical resolution 
for band lA iheech et al. 2002|). 
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Fig. A. 5. Comparison between the AOT06 observation 
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Fig. A.3. Band 2 for Sirius (Al V) and 7 Dra (K5 III) 

are displayed. The ISO-SWS spectrum is plotted in black, 
the synthetic spectrum in red. Problems with the RSRF 
are indicated by an arrow. 
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Fig. A. 6. Top: Comparison between the AOTOl speed- 
4 observation of 7 Dra rebinned at 1) a resolving power 
of 1500 (black) and 2) a resolving power of 1300 (red). 
Bottom: Comparison between the ISO-SWS AOTOl ob- 
servation in band lA (black) and its synthetic spectrum 
(red) at a resolving power of 1300. 
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Fig. A. 7. Summer (red) and winter (green) FTS-KP 
spectra of a Boo at a resolving power A/AA of 60000. 
They are compared with the synthetic spectrum (black) 
for a Boo with stellar parameters Toff ~ 4350 K, log g 
= 1.50, M = 0.75Mo, [Fe/H] = -0.50, 6 = 1.7km s-\ 
12C/13C = 7, £(C) = 7.96, e(N) = 7.61, £(0) = 8.68, 
e(Mg) = 7.33 and Ft = 3km s^^. The OH 1-0 lines are 
predicted as too weak. 
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1. Introduction 

For the astronomical community analysing ISO-SWS 
data (Infrared Space Observatory, Short- Wavelength 
Spectrometer, de Graauw et al. 1996), a first point to 
assess when judging and qualifying their data concerns 
the flux calibration accuracy. Since the calibration process 
is not straightforward, knowledge on the full calibration 
process and on the still remaining calibration problems is 
crucial when processing the data. 

One way to detect calibration problems is by compar- 
ing observed data with theoretical predictions of a whole 
sample of standard calibration sources. But, as explained 
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in Decin ct al. (2000) (hereafter referred to as Paper I) 
a full exploitation of the ISO-SWS data may only re- 
sult from an iterative process in which both new theo- 
retical developments on the computation of stellar spec- 
tra — based on the MARCS and Turbospectrum code 
(Gustafsson et al. 1975; Plez et al. 1992, 1993), version 
May 1998 — and more accurate instrumental calibration 
arc involved. 

Precisely because this research entails an iterative pro- 
cess, one has to be extremely careful not to confuse techni- 
cal detector problems with astrophysical issues. Therefore, 
the analysis in its entirety encloses several steps. Some 
steps have already been demonstrated in the case of a 
Tau in Paper I. They will be summarised in Sect. 2. Other 
points will be introduced in Sect. 2 and will be elaborated 
on in the first sections of this article (Sect. 3-5). Having 
described the method of analysis, the general discrepancies 
between observed and synthetic spectra are subjected to a 
careful scrutiny in order to elucidate their origin. At this 
point, a distinction can be made between discrepancies 
typically for warm stars and those typical for cool stars. 
For this research, warm stars arc defined as being hotter 
than the Sun (Teff,© = 5770 K) and their infrared spec- 
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tra are mainly dominated by atomic lines, while molecular 
lines are characteristic of cool star spectra. A description 
on the general trends in the discrepancies for warm and 
cool stars will be made in this paper, while each star of the 
sample will be discussed individually in two forthcoming 
papers in which also an overview of other published stellar 
parameters will be given. 

As stated in Paper I. the detailed spectroscopic analy- 
sis of the ISO-SWS data has till now been restricted to the 
wavelength region from 2.38 to 12 ^m. So, if not specified, 
the wavelength range under research is limited to band 1 
(2.38 - 4.08 Afm) and band 2 (4.08 - 12.00 /im). Band 3 
(12.00 - 29.00 i^Lui) will be elaborated on by Van Malderen 
(Van Malderen ct al., 2001, in prep.). 

This paper is organised as follows: in Sect. 2 the gen- 
eral method of analysis is summarised. The sample of ISO- 
SWS observations is described in Sect. 3, while the data 
reduction procedure is discussed in Sect. 4. The obser- 
vations of two independent instruments are introduced in 
Sect. 5. In Sect. 6, the results are elaborated on. In the last 
section, Sect. 7, the impact on the calibration of ISO-SWS 
and on the theoretical modelling is given. 

The appendix of this article is published electronically. 
Most of the grey-scale plots in the article are printed in 
colour in the appendix, in order to better distinguish the 
different spectra. 

2. Description of the strategy 

Since this research includes an iterative process, one has 
to be very careful not to introduce (and so to propagate) 
any errors. The strategy in its entirety therefore encom- 
passes a number of steps, including (1.) a spectral cov- 
erage of standard infrared sources from AO to M8, (2.) 
a homogeneous data reduction, (3.) a detailed Hterature 
study, (4.) a detailed knowledge of the impact of the vari- 
ous parameters on the spectral signature, (5.) a statistical 
method to test the goodncss-of-fit (Kolmogorov-Smirnov 
test) and (6.) high-resolution observations with two inde- 
pendent instruments. Some points (in particular, points 4 
and 5) have already been demonstrated in the case of a 
Tau in Paper I. Points 1, 2 and 6 will be elaborated on in 
Sect. 3, Sect. 4 and Sect. 5 respectively. 

In its totality, the general method of analysis — based 
on these 6 points — may be summarised as follows: 
a large set of standard stars (AO - M8) has been ob- 
served with ISO-SWS (Sect. 3). The observational data 
first have been subjected to a homogeneous data-reduction 
procedure (Sect. 4). Thereafter, the carefully reduced ISO- 
SWS data of one warm and one cool star were compared 
with the observational data of two independent instru- 
ments (FTS-KP and FTS-ATMOS, see Sect. 5). This step 
is very crucial, since this is the only secure and decisive 
way to point out calibration problems with the detectors of 
ISO-SWS. The complete observational data-set, covering 
a broad parameter space, was then compared with theo- 
retical predictions. By knowing already some problematic 
points in the cahbration of ISO-SWS, these comparisons 



led both to a refinement of our knowledge on the cali- 
bration problems and to a determination of theoretical 
modelling problems (Sect. 6). 

The knowledge on the relative importance of the 
different molecules ^ — displaying their characteristic 
absorption pattern somewhere in the broad ISO-SWS 
wavelength-range — and on the impact of the various stel- 
lar parameters on the infrared spectrum enabled us also 
to determine the fundamental stellar parameters for the 
cool giants in our sample (Paper I). Due to severe cali- 
bration problems with the band-2 data (see Sect. 4), only 
band-1 data were used for this part of the process. Once 
a high-level of agreement between observed and synthetic 
data was reached, a statistical test was needed to objec- 
tively judge on the different synthetic spectra. A choice 
was made for the Kolmogorov-Smirnov test (Paper I). 
This statistical test globally checks the goodness-of-fit of 
the observed and synthetic spectra by computing a devia- 
tion estimating parameter j3 (see Eq. (5) in Paper I). The 
lower the /?- value, the better the accordance between the 
observed data and the synthetic spectrum. 

Using this method, the effective temperature, grav- 
ity, metallicity, microturbulent velocity together with the 
abundance of C, N and O and the ^^C/^'^C-ratio were es- 
timated for the cool stars. From the energy distribution 
of the synthetic spectrum between 2.38 and 4.08 /xm and 
the absolute flux-values in this wavelength range of the 
ISO-SWS spectrum, the angular diameter was deduced. 
We therefore have minimised the residual sum of squares 

n 

E(/«-5«)', (1) 

i=l 

with f{i) and g{i) representing respectively an observa- 
tional and a synthetic data point at the zth wavelength 
point. 

The error bars on the atmospheric parameters were es- 
timated from 1. the intrinsic uncertainty on the synthetic 
spectrum (i.e. the possibility to distinguish different syn- 
thetic spectra at a specific resolution, i.e. there should be 
a significant difference in /3-values) which is thus depen- 
dent on both the resolving power of the observation and 
the specific values of the fundamental parameters, 2. the 
uncertainty on the ISO-SWS spectrum which is directly 
related to the S/N of the ISO-SWS observation, 3. the 
value of the /3-parameters in the Kolmogorov-Smirnov test 
and 4. the still remaining discrepancies between observed 
and synthetic spectra. 

^ As for atoms, laboratory data for molecules are generally 
more accurate than computed data, but much too sparse to 
serve as the sole source for the computation of photospheric 
models or to represent the myriads of weak lines in a cool 
star spectrum. An overview of existing databases can be found 
in AUard et al. (1992); Grevesse & Sauval (1992); Gustafsson 
& Jorgensen (1994); j0rgensen et al. (1992); Kurucz (1992); 
Decin (2000). Some of these, useful for our research, are de- 
scribed in Decin (2000). A copy of this description can be re- 
trieved from http;//www.ster.kuleuven.ac.be/'leen. The line 
lists used for this research are enumerated in Paper I. 
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It should be noted that an error on the effective 
temperature introduces an error on the angular diame- 
ter. The IR flux of the cool giants does not follow the 
Rayleigh- Jeans law for a black-body, and we can write 

oc Teff« • ej. Thus, with oc ^^^3'^' ■ Teff"«/' and 
o'<jF> = 0.10J-, one obtains 

with cr?,^ being the intrinsic uncertainty on the angular di- 
ameter (i.e. the amount by which one may change the an- 
gular diameter without significant difference in the resid- 
ual sum of squares). The uncertainty on the angular diam- 
eter (Eq. (2)) is mainly determined by the uncertainty in 
the absolute flux (the first term in Eq. (2)), resulting in al- 
most the same percentage errors in Bd - Using the ISO-SWS 
observational data of these cool giants, we could demon- 
strate that q « 1.3 ± 0.1 for 3600 K < Tcff< 4400 K 
(Decin 2000). With q « 1.3, cr<T„a> = 70K and Toff = 
3850 K, the uncertainty increases with a factor 1.008 com- 
pared to the q ~ I situation. 

From the angular diameter and the parallax measure- 
ments (mas) from Hipparcos (with an exception being a 
Cen A, for which a more accurate parallax by Pourbaix 
et al. 1999, is available), the stellar radius was derived. 
This radius, together with the gravity — determined 
from the ISO-SWS spectrum — then yielded the gravity- 
inferred mass. From the radius and the effective tempera- 
ture, the stellar luminosity could be extracted. 

This method of analysis could however not be applied 
to the warm stars of the sample. Absorption by atoms de- 
termines the spectrum of these stars. It turned out to be 
unfeasible to determine the effective temperature, gravity, 
microturbulence, metallicity and abundances of the chem- 
ical elements from the ISO-SWS spectra of these warm 
stars, due to 

— problems with atomic oscillator strengths in the in- 
frared (see Sect. 6.1), 

— the small dependence of the IR continuum on the fun- 
damental parameters (when changed within their un- 
certainty) , 

— the small dependence of the observed atomic line 
strength on the fundamental parameters, and 

— the absence of molecules, which are each of them specif- 
ically dependent on the various stellar parameters. 

Therefore, good-quality published stellar parameters 
were used to compute the theoretical model and corre- 
sponding synthetic spectrum. The angular diameter was 
deduced directly from the ISO-SWS spectrum, which then 
yielded — in conjunction with the assumed parallax, grav- 
ity and effective temperature — the stellar radius, mass 
and luminosity. 

3. Description of the sample 

It seems important to cover a broad parameter space in or- 
der to be able to distinguish between calibration problems 



and problems related to the computation of a theoretical 
model and/or to the generation of a synthetic spectrum. 
Therefore, stellar standard candles spanning the spectral 
types AO - MS were observed. The observations were ob- 
tained in the context of two ISO-SWS proposals. Some 
calibration data were also provided by the SIDT (SWS 
Instrument Dedicated Team) in the framework of a quick- 
feedback refining of the model SEDs used for the SWS 
calibration. In total, full-resolution SWS scans of 20 stan- 
dard stars, covering the full AO - M8 spectral range were 
obtained. Stars with an earlier spectral type were included 
in the proposals of S. Price. An overview of the objects, ob- 
serving dates and integration times (tint) can be found in 
Tables 1-2. Stars indicated by a '•' are stars which have 
been scrutinised more carefully than stars indicated by a 
'o'. We first focused on a Bootes (K2 IIIp), because of its 
well-known stellar parameters and the high quality of the 
ISO-SWS data for this object (Decin et al. 1997). From 
then on, we have gone, step by step, towards both higher 
and lower temperatures. Why the two warmest stars in 
our sample (being Vega and Sirius) have been considered, 
but have not been studied into all detail, will be explained 
in Sect. 6.1. The other 'o '-stars in Table 2 are stars cooler 
than an M2 giant — i.e. cooler than ~ 3500 K. Calibration 
problems with 7 Cru, variability (Monnier et al. 1998), 
the possible presence of a circumstellar envelope, stellar 
winds or a warm molecular envelope above the photo- 
sphere (Tsuji et al. 1997) made the use of hydrostatic 
models for these stars inappropriate and have led to the 
decision to postpone the modelling of the coolest stars in 
the sample. 

4. Data reduction 

In order to reveal calibration problems, the ISO-SWS data 
have to be reduced in a homogeneous way. For all the 
stars in our sample, at least one AOTOl observation^ 
(AOT = Astronomical Observation Template; AOTOl = a 
single up-down scan for each aperture with four possible 
scan speeds at degraded resolution) is available, some stars 
have also been observed using the AOT06 mode ( = long 
up-down scan at full instrumental resolution). Since these 
AOTOl observations form a complete and consistent set, 
they were used as the basis for the research. In order to 
check potential calibration problems, the AOT06 data are 
used. The scanner speed of the highest-quality AOTOl 
observations was 3 or 4, resulting in a resolving power ~ 
870 or ~ 1500, respectively (Leech et al. 2002). The ap- 
propriate resolving power of each sub-band was taken to 
be the most conservative theoretical resolving power as 
determined by Lorcnte in Leech et al. (2002), with the 
exception being band lA for which this value has been 

^ Each observation is determined uniquely by its observation 
number (8 digits), in which the first three digits represent the 
revolution number. The observing data can be calculated from 
the revolution number which is the number of days after 17 
November f995. 
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Table 1. The logbook of the ISO-SWS observations used for 
this study. The observing date can be calculated from the rev- 
olution number which is the number of days after 17 November 
1995. Sub-bands of bands 1 and 2 observed completely or par- 
tially during an AOT06 observation are given underneath the 
table. The wavelength ranges for the different sub-bands can 
be found in Table 3. 



name 


HD 


HR 


HIC 


Spectral 


AOT mode 




revo- 










TvDe 


(speed) 


['"'] 


lution 


a Lyr 


172167 


7001 


91262 


AOV 


AOTOl (3) 


3462 


178 


Vega 










AOT06 


5642 


650*^ 












AOT06 


4354 


678*^ 


a CMa 


48915 


2491 


32349 


AlV 


AOTOl (4) 


6538 


689 


Sirius 










AOTOl (1) 


1140 


868 


• (3 Leo 


102647 


4534 


57632 


A3Vv 


AOTOl (3) 


3462 


189 


Denebola 










AOTOl (1) 


1096 


040 


• a Car 


'±oo'±o 


2326 




Foil 




uOoo 


729 


Canopus 














• a Ceil A 


128620 


5459 


71683 


G2V 


AOTOl (4) 


6538 


607 












AOTOl (1) 


1140 


294 


• 5 Dra 


180711 


7310 


94376 


G9III 


AOTl (4) 


6538 


206 












AOTOl (4) 


6528 


072 


• ^ Dra 


163588 


6688 


87585 


K2III 


AOTOl (3) 


3454 


314 












AOTOl (1) 


1044 


068 


• a Boo 


124897 


5340 


69673 


K2IIIp 


AOTOl (4) 


6538 


452 


Arcturus 










AOTOl (4) 


6528 


071 












AOTOl (1) 


1140 


275 












AOTOl (1) 


1094 


056 












AOT06 


3904 


583*^ 












AOT06 


4720 


601*" 












AOT06 


4510 


608*^^ 


• Q Tuc 


211416 


8502 


110130 


K3III 


AOTOl (4) 


6539 


866 


• (3 UMi 


131873 


5563 


72607 


K4III 


AOTOl (4) 


6546 


182 


Kochab 










AOTOl (2) 


1816 


079 



*i: lA, IB, ID, IE, 2A, 2B, 2C; lA, ID, IE, 2A; lA, 
IB, ID, 2A, 2C; lA, IB, ID, IE, 2A, 2B, 2C; lA, IB, 
ID, IE, 2A, 2B, 2C 



changed from 1500 to 1300, as will be discussed in Sect. 
6.2. 

The ISO-SWS data were processed to a calibrated 
spectrum by using the same procedure as described in 
Paper I using the calibration files available in OLP6.0. 

The band 2 (Si:Ga) detectors used in SWS 'remember' 
their previous illumination history. Going from low to high 
illumination, or vice versa, results in detectors asymptot- 
ically reaching their new output value. These are referred 
to as memory effects or transients. For sources with fluxes 
greater than about 100 Jy, memory effects cause the up 
and down scans in the SPD ( = Standard Processed Data) 
to differ in response by up to 20% in band 2. Since an 



Table 2. Logbook of the ISO-SWS observations: continuation 
of Table 1. 



name 


HD 


HR 


HIC 


opectrai 


± nioQC 


J mt 












Type 


(speed) 


[sec] 


lution 


• 7 Dra 


164058 


6705 


87833 


K5III 


AOTOl (4) 


6538 


377 


Etamin 










AOTOl (4) 


6542 


040 












AOTOl (2) 


1912 


811 












AOTOl (1) 


1140 


496 












AOTOl (1) 


1062 


126 












AOT06 


4568 


501*"^ 












AOT06 


5676 


559*^ 


• Q Tau 


29139 


1457 


21421 


K5III 


AOTOl (4) 


6538 


636 


Aide bar an 










AOT06*** 


3700 


681 


• H See 
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adapted version of the Fouks-Schubert model to correct 
for these memory effects in band 2 was still in develop- 
ment (Leech et al. 2002), this method could not be applied 
during our reduction procedure. Instead, we have used the 
down-scan data of our observation as a reference to do a 
correction of the flux level of the flrst scan (up-scan). This 
is justifled since the memory effects appear to be less se- 
vere in the down-scan measurements, suggesting a more 
stabilised response to the flux level for the down-scan data. 

Also the band 2 dark current subtraction is closely 
tied to the band 2 memory effect correction. The mem- 
ory effect for Si:Ga detectors as described by the Fouks- 
Schubert model is an additive effect. As such, its proper 
correction will take place during the dark current subtrac- 
tion. Since this correction tool was still not available, all 
dark currents were checked individually. When a dark cur- 
rent was corrupted too much by memory effects, its value 
was replaced by the value of a preceding or following dark- 
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current not being affected. In this way, a small error can 
occur, which is, however, negligible due to the high flux 
level of our stellar sources. 

Table 3. Factors used to multiply the sub-bands for the 
highest- quality AOTOl observation (denoted by its revolution 
number) of the selected stars in our sample are given. For each 
band the wavelength range Xt - Xe (in /xm) and the (con- 
stant) resolving power for a speed-4 AOTOl observation are 
also noted. For a speed-3 AOTOl observation one has to mul- 
tiply the average resolving power by ~ 0.58. In the last two 
columns the pointing offset as estimated from known imper- 
fections in the satellite attitude control system is given in arc- 
sec, where the y-axis denotes the cross-dispersion direction for 
SWS and the z-axis the dispersion direction. 



lA IB ID IE 2A 2B 2C rev. dy dz 



A/ A A 1300 1200 1500 1000 1200 800 800 
Ai, [pm] 2.38 2.60 3.02 3.52 4.08 5.30 7.00 
Ae [fim] 2.60 3.02 3.52 4.08 5.30 7.00 12.00 



a Lyr 1.06 1.06 1.00 1.00 1.00 0.97 +12Jy 178 -0.608 -1.179 



a CMa 1.00 1.00 1.00 0.995 1.12 1.23 1.00 689 0.034 0.003 



/3 Leo 0.99 0.99 1.00 1.00 1.16 1.27 +3.5Jy 189 0.478 0.556 



a Car 0.97 0.98 1.00 1.00 0.98 1.10 0.91 729 0.024 0.072 



a Ccn A 1.01 1.02 1.00 1.01 0.985 1.06 0.91 607 0.000 0.000 



S Dra 0.97 0.98 1.00 1.015 1.03 1.02 1.10 206 -0.422 1.480 



5 Dra 0.99 0.99 1.00 0.99 1.12 1.15 1.05 314 1.286-0.282 



a Boo 0.995 1.01 1.00 1.005 0.95 1.05 1.00 452 0.000 0.000 



a Tuc 1.005 1.02 1.00 1.01 1.05 1.00 1.00 866 0.000 0.000 



/3 UMi 1.00 1.015 1.00 1.01 0.91 0.885 1.00 182 -1.062 0.045 



7 Dra 0.995 1.005 1.00 1.005 0.935 0.98 0.91 377-0.304 0.181 



a Tau 1.00 1.01 1.00 1.00 1.00 1.045 1.00 636 0.000 0.000 



H Sco 1.00 1.015 1.00 1.00 1.13 1.05 1.15 847 0.000 0.000 



And 1.00 1.00 1.00 1.005 1.00 1.10 0.95 795 0.000 0.000 



a Cot 0.985 1.00 1.00 1.01 0.935 1.03 0.91 797 0.000 0.000 



/3 Peg 1.00 1.015 1.00 1.005 0.935 1.03 0.935 551 0.017 0.095 



Since the different sub-band spectra can show jumps 
in flux at the band-edges, several sub-bands had to be 
multiplied by a small factor to construct a smooth spec- 
trum. Three causes for the observed shift factors between 



different sub-bands of an observation and between differ- 
ent observations of a given stellar source can be reported: 
1. pointing errors, 2. problems with the RSRF correction, 
and 3. a problematic dark current subtraction, from which 
the pointing errors are believed to have the largest impact. 
The pointing errors as well as the RSRF correction causes 
a decrease in flux by a gain factor, while the dark current 
subtraction can lower the flux level by an offset. As the 
effects of the pointing errors are estimated to have the 
biggest effect, and since the stars in our sample have a 
high flux level so that the dark current subtraction only 
plays a marginal role, the individual sub-bands were multi- 
plied with a factor — rather than shifted with an offset — 
in order to obtain a smooth spectrum. These factors (see 
Table 3) were determined by using the overlap regions of 
the different sub-bands and by studying the other SWS 
observations. The band-ID data were taken as reference 
data, due to the absence of strong molecular absorption 
in this wavelength range which may cause a higher stan- 
dard deviation in the bins obtained when rebinning the 
oversampled spectrum, and — most importantly — due 
to the low systematic errors in this band, caused by e.g. 
errors in the curve of the RSRF, detector noise, uncertain- 
ties in the conversion factors from fiV/s to Jy, ... (Leech 
et al. 2002). Using the total absolute uncertainty values 
— which have accumulated factors from each of the cali- 
bration steps plus estimated contributions from processes 
which were unprobed or uncorrected — as given in Table 
5.3 in Leech et al. (2002), the estimated la uncertainty 
on these factors is 10%. As is clearly visible from Table 
3, these factors do not show any trend with spectral type 
or flux- level. This is displayed in Fig. 1, where the band- 
border ratios between lA-lB, IB-ID and ID-IE are plot- 
ted in function of the flux at 2.60, 3.02 and 3.52 //m respec- 
tively. For this plot, all the observations of the cool stars in 
our sample, discussed in the Appendix of Paper IV of this 
series, are used. In band 1, the band-border ratios of lA- 
IB and ID-IE are from bands within the same aperture. 
Going from band IB to band ID, the aperture changes. 
Satellite mispointings can have a pernicious impact on this 
band-border ratio: the mean deviation of the band-border 
ratios w.r.t. 1 is signiflcantly larger for IB-ID {— 0.015) 
than for lA-lB and ID-IE (being respectively 0.009 and 
0.005). Due to the problems with memory effects in band 
2 (4.08 - 12 /im), the factors of each sub-band of band 
2 were determined by use of the corresponding spectral 
data of Cohen (Cohen et al. 1992, 1995, 1996; Witteborn 
et al. 1999): for Vega and Sirius Cohen has constructed 
a calibrated model spectrum; a composite spectrum (i.e. 
various observed spectra have been spliced to each other 
using photometric data) is available for a Cen A, a Boo, 
7 Dra, a Tau, /3 And, a Cet, and (3 Peg; a template spec- 
trum (i.e. a spectrum made by using photometric data of 
the star itself and the shape of a 'template' star) is built 
for 6 Dra (template: (3 Gem: KO III), ^ Dra (template: a 
Boo: K2 IIIp), a Tuc (template: a Hya: K3 ILIII) and 
H Sco (template: a Tau: K5 III). When no template was 
available (for (3 Leo, a Car and /3 UMi), the synthetic spec- 
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trum showing the best agreement with the band-1 data 
was used as reference. This does not imply that we are 
trapped in a circular argument, since the stellar parame- 
ters for the synthetic spectrum were determined from the 
band-1 data only. Moreover, the maximum difference in 
the correction factors for band 2 obtained when using the 
synthetic spectra instead of a Cohen template for the 13 
stars common in the sample is 7 %, which is well within the 
photometric absolute flux uncertainties claimed by Leech 
et al. (2002). Note that all shift factors are in within the 
AOTOl band border ratios as derived in Fig. 5.33 and Fig. 
5.34 in Leech et al. (2002). Using the overlap regions in 
band 2 can have quite a big effect on the final composed 
spectrum: focussing on f3 UMi, we note that by using these 
overlap regions band 2A (and consequently bands 2B and 
2C) should be shifted downwards by a factor 1.04; in order 
then to match the shifted band 2B and band 2C, band 2C 
should be once more shifted downwards by a factor 1.12. 
In general, the error in the absolute flux could increase to 
~ 20 % at the end of band 2C when this method would be 
used. 

For a more elaborate discussion on the SWS error bud- 
get, we would hkc to refer to Leech et al. (2002). 
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Fig. 1. Flux ratios between the different sub-bands of band 1 
at the wavelengths of overlap. The mean deviation w.r.t. 1 is 
given for the different band borders at the right corner of the 
figure. 



5. Verification of tlie results witli FTS-KP and 
FTS-ATMOS spectra 

To test our findings — indicating either a calibration 
problem or a theoretical modelling problem — with 
data taken with an independent instrument, a high- 
resolution observation of both one warm and one cool 
star were included. The high-resolution Fourier Transform 
Spectrometer (FTS) spectrum of a Boo (Hinkle et al. 
1995b) and the Atmospheric Trace Molecule Spectroscopy 



(ATMOS) spectrum of the Sun (Farmer & Norton 1989; 
Geller 1992) are used as external control to the process. 

The Arcturus observations were obtained with the 
FTS at the Kitt Peak (KP) 4m MayaU telescope mainly in 
1993 and 1994 and are described in Hinkle et al. (1995a, b). 
The entirety of the 1 to 5 fim Arcturus spectrum de- 
tectable from the ground was observed twice at opposite 
heliocentric velocities, allowing the removal of many tel- 
luric lines. We refer to these two spectra as the winter and 
summer observations. The data were obtained at a resolv- 
ing power (A/AA) near 100000. Spectra of the Earth's 
atmosphere, derived from high resolution solar spectra, 
have been ratioed to the Arcturus spectra to largely re- 
move the telluric lines. Most telluric lines less than 30% 
deep are cleanly removed but as the optical depth of the 
telluric lines increases more information is lost and lines 
over 50 % deep can not be removed. Gaps in the plots ap- 
pear at these spectral regions. A typical example of the 
Arcturus spectra is shown in Fig. 16 on page 15. In the 
upper panel of Fig. 16, some CO lines of the first overtone 
arc plotted. In the lower panel of Fig. 16, one can clearly 
distinguish atomic and OH lines. 

During the period from November 3 to November 
14, 1994, the Atmospheric Trace Molecule Spectroscopy 
(ATMOS) experiment operated as part of the ATLAS-3 
payload of the shuttle Atlantis (Gunson et al. 1996). The 
principal purpose of this experiment was to study the dis- 
tribution of the atmosphere's trace molecular components. 
The instrument, a modified Michelson interferometer, cov- 
ering the frequency range from 625 to 4800 cm~^ at a spec- 
tral resolution of 0.01 cm~^, also recorded high-resolution 
infrared spectra of the Sun. A small part of the inten- 
sity spectrum is shown in Fig. 2. The Holweger-Miiller 
model (Holweger & Miiller 1974) was used as input for 
the computation of the synthetic spectrum of the Sun. A 
microturbulent velocity of 1 km s~^ was assumed. 

In the case of a Boo, a high-resolution synthetic spec- 
trum was first generated with the stellar parameters of 
Kja^rgaard et al. (1982) as input parameters (Teg = 
4350 K, log g = 1.80, [Fe/H] ^ -0.50, e(C) = 7.89, e(N) = 
7.61, e(0) = 8.68). For the microturbulence a velocity of 
1.7 km s^^ was assumed as being the mean for red giants 
(Gustafsson et al. 1974). Since both Mackle et al. (1975) 
and Bell et al. (1985) found a Si and Mg abundance de- 
pleted with respect to the solar values, e(Mg) was taken to 
be 7.33 and £(Si) to be 7.20. For the anisotropic macrotur- 
bulence Ft, a radial-tangential profile was assumed (Gray 
1975) with FWHM of 3 km s'^. This synthetic spectrum 
was then convolved with the beam-profile function (i.e. 
a sine- function) . Both observed FTS-KP and synthetic 
spectrum were rcbinncd to a resolving power of 60000, 
since the resolving power of the FTS-KP was not con- 
stant over the whole wavelength range. Using this set of 
input-parameters, the high-excitation CO (Aw = 1 and 
Av = 2) and SiO {Av = 2) lines were predicted as be- 
ing too weak. This was confirmed from both the FTS-KP 
and the ISO-SWS spectrum. Especially from the ISO-SWS 
spectrum — and more specifically from the slope in bands 
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parameters will be discussed and confronted with other 
published stellar parameters. 
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Fig. 2. FTS-ATMOS intensity spectrum of the Sun (black) 
at a resolving power A/ A A of 60000. The synthetic spectrum 
(grey) has been computed using the Holweger-Miiller model 
(Holweger & Miiller 1974) with = 1km s"^ 



IB and ID — it became clear that the gravity should be 
lowered. Using then the FTS-KP-spectra, the input pa- 
rameters were improved until an optimal fit was obtained 
(see Fig. 16 on page 15). This resulted in the following 
parameters: Teff = 4350 K, log g = 1.50, [Fe/H] = -0.50, 
= 1.7km s-i, 12C/13C = 7, e(C) = 7.96, e(N) = 7.61, 
e(0) = 8.68, £(Mg) = 7.33 and e(Si) = 7.20. 

Discrepancies appearing from the confrontation of the 
synthetic spectrum with both the FTS-KP and the ISO- 
SWS spectra, may be clearly attributed to problems in 
constructing the model or generating the synthetic spec- 
trum. Other discrepancies in the ISO-SWS versus syn- 
thetic spectrum are the ones caused by calibration prob- 
lems. 

For the warm stars in the sample, an analogous com- 
parison was performed with the help of the FTS-ATMOS 
intensity spectrum of the Sun. Discrepancies between the 
ISO-SWS and synthetic spectra of a Gen A (G2V) were 
compared with the discrepancies found in the comparison 
between the FTS-ATMOS spectrum of the Sun and its 
synthetic spectrum (Fig. 2). 

Inspecting Fig. 2 and Fig. 16, some discrepancies be- 
come immediately visible, e.g. the CO 2-0 and 3-1 are pre- 
dicted too weak in both figures. The reason for this and 
other discrepancies will be explained in the next section. 

6. Results 

Computing synthetic spectra is one step, distilling useful 
information from it is a second — and far more difficult 
— one. Fundamental stellar parameters for this sample of 
bright stars are a first direct result which can be deduced 
from this comparison between ISO-SWS data and syn- 
thetic spectra. In papers III and IV of this series, these 
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Fig. 3. Comparison between the ISO-SWS data of a Cen A 
(black) and the synthetic spectrum (grey) with stellar parame- 
ters Tcft = 5830 K, log g = 4.35, M = 1.3 M©, [Fe/H] = 0.25, Ct 
= 1.0kms-\ i^c/i^C = 89, e(C) = 8.74, e(N) = 8.26, e(0) = 
9.13 and 9d ~ 8.80 mas. Some of the most prominent discrep- 
ancies between these two spectra are indicated by an arrow. 
A coloured version of this plot is available in the Appendix as 
Fig. ??. 



A typical example of both a warm and cool star is 
given in Fig. 3 and Fig. 4 respectively. Different types of 
discrepancies do emerge. The size of the discrepancies be- 
tween ISO-SWS observations and theoretical predictions 
varies a lot. Both for the warm and for the cool stars we 
see a general good agreement in shape between the ISO- 
SWS and theoretical data in band 1 , with however (local) 
error peaks up to ~ 8 %. An exception is f3 Peg for which 
the shape is wrong by ~ 6 % and local error peaks may go 
up to ~ 15 %. The agreement between observational and 
synthetic data is worse in band 2: a general mismatch by 
up to ~ 15 % may occur. 

By scrutinising carefully the various discrepancies be- 
tween the ISO-SWS data and the synthetic spectra of the 
standard stars in our sample, the origin of the different 
discrepancies was elucidated. First of all, a description on 
the general trends in discrepancies for the warm stars will 
be made, after which the cool stars will be discussed. 

6.1. Warm stars: AO - G2 

1. Especially when concentrating on a Cen A (G2 V), 
one notifies quite a few spectral features which ap- 
pear in the ISO-SWS spectrum, but are absent in the 
synthetic spectrum. By comparing the spectra of the 
other warm stars {f3 Leo, a Car, a Lyr and a CMa) 
with each other, corresponding spectral features can 
be recognised in their ISO-SWS spectra, although this 
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Fig. 4. Comparison between band 1 and band 2 of the ISO- 
SWS data of 7 Dra (black) and the synthetic spectrum (grey) 
with stellar parameters Toff — 3960 K, log g — 1.30, M — 
1.7Mq, [Fe/H] = 0.00, Ct = 2.0km s-\ ^^c/i^c = 10, e(C) = 
8.15, £(N) = 8.26, e(0) = 8.93 and Od = 9.98 mas. A coloured 
version of this plot is available in the Appendix as Fig. ??. 



is somewhat more difficult for (3 Leo due to the lower 
resolving power and lower signal-to-noise ratio. Also 
for the cooler stars in the sample, these spectral fea- 
tures are (weakly) present. Some of the most promi- 
nent ones are indicated by arrows in the spectrum of 
a Cen A in Fig. 3. Molecular absorption can not be a 
possible cause of/ contribution to e.g. the 3.9 /im fea- 
ture, since this feature is seen for the warm stars in the 
sample, where no molecular absorption occurs at these 
wavelength ranges. The solar FTS-ATMOS spectrum 
proved to be extremely useful for the determination of 
the origin of these features. All spectral features, indi- 
cated by an arrow in Fig. 3, turned out to be caused by 
— strong — atomic lines (Mg, Si, Fc, Al, C, ...). This is 
illustrated for the wavelength range from 3.85 to 3.92 
/im in Fig. 5. In panel (a) of Fig. 5, the high-resolution 
FTS-ATMOS spectrum of the Sun is compared with its 
synthetic spectrum in the wavelength range from 3.853 
to 3.917 yum. The strongest lines are identified by using 
the line list of Geller (1992). At the ISO resolving power 
of 1000 in band IE, these atomic lines are reduced to 
the features indicated in panel (b). The same features 
can be recognised in the ISO-SWS spectrum of a Cen 
A in panel (c) . It is clear that these atomic features are 
not well calculated for the synthetic spectra of the Sun 
and a Cen A. 
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Fig. 5. In panel (a) the high-resolution FTS-ATMOS spectrum 
of the Sun is compared with its synthetic spectrum based on 
the Holweger-Miiller model (1974) and computed by using the 
atomic line list of Hirata & Horaguchi (1995). Panel (b) shows 
the same comparison as panel (a) , but now at a resolving power 
of 1000 for the wavelength range going from 3.52 - 4.08 fj.m 
(band IE). The data of band IE of the ISO-SWS spectrum of 
a Cen A (at a resolving power of 1000) are plotted in panel 
(c). 



For Fig. 5, the atomic line list of Hirata & Horaguchi 
(1995) was used to generate the synthetic spectrum. 
This line list has as starting files the compilation by 
Kurucz & Peytremann (1975) and Kurucz (1989) who 
lists the semi-empirical (^/-values for many ions. Energy 
levels were adopted from recent compilations (Sugar 
& Corliss 1985) or individual works. The line list of 
Kelly (1983) was merged into the file. The gf-values 
are taken from several compilations (e.g. Fuhr et al. 
1988a,b; Reader et al. 1980; Wiese et al. 1966, 1969; 
Morton 1991, 1992). Published results of the Opacity 
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Project (Seaton 1995) were also included. Comparing 
the atomic line list of Hirata & Horaguchi (1995) with 
the identifications as given by GcUcr (1992) made clear 
that quite some lines are misidentified or not included 
in the IR atomic line list of Hirata & Horaguchi (1995). 
The usage of VALD (Vienna Atomic Line Database: 
Piskunov ct al. 1995; Ryabchikova et al. 1997; Kupka 
et al. 1999)3 and the line list of van Hoof (1998) did not 
solve the problem. Obviously, the oscillator strengths 
of the atomic lines in the infrared arc not known suffi- 
ciently well. 

In order to test this hypothesis, Sauval (2002, priv. 
comm.) has constructed a new atomic line list by 
deducing new oscillator strengths from the high- 
resolution ATMOS spectra of the Sun (625 - 4800 
cm~^). A preliminary comparison using this new 
linelist is given in Fig. 6 (which should be compared 
with Fig. 5). The contents, nature, limitations, uncer- 
tainties, ... of this new line list will be discussed in 
a forthcoming paper. But it is already obvious from a 
confrontation between Fig. 5 and Fig. 6 that these new 
oscillator strengths from Sauval are more accurate. 
It is plausible that the same reason, i.e. wrong and miss- 
ing oscillator strengths of atomic lines in the infrared, 
together with noise, is the origin of the observed dis- 
crepancies between the ISO-SWS and synthetic spectra 
for the other warm stars in the sample, because 

— incomplete atomic data already caused a problem 
for the synthetic spectrum of a Cen A, and 

— the stars concerned are hotter than a Cen A, so 
atoms (maybe of higher excitation or ionisation) 
will determine even more the spectral signature of 
these stars. 

The lack of reliable atomic data rendered the deter- 
mination of the fundamental stellar parameters for 
the warm stars from the ISO-SWS data impossible. 
Another consequence concerns the continuum, which 
was very difficult to determine. Therefore, the uncer- 
tainty on the angular diameter of these warm stars is 
more pronounced. 
2. Secondly, the hydrogen lines are conspicuous. For ex- 
ample, the synthetic hydrogen Pfund lines are almost 
always predicted as being too strong for main-sequence 
stars, while they are predicted as being too weak for 
the supergiant a Car. This indicates a problem with 
the generation of the synthetic hydrogen lines, which 
is corroborated when the high-resolution FTS- ATMOS 
spectrum of the Sun is compared with its synthetic 
spectrum (sec Fig. 7). In the TurboSpcctrum program 
(Plez et al. 1992, 1993) the hydrogen line-profile cal- 
culation, adopted from the SYNTHE code of Kurucz, 
includes the Stark, the Doppler, the van dcr Waals and 
the resonance broadening. For the first four lines in cv- 



^ No more recent "near IR" data for neutral metal lines have 
been added to VALD, so no improvements for this particular 
wavelength range could be expected over the Kurucz (1992) 
line lists. 
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Fig. 6. In panel (a) the high-resolution ATMOS spectrum of 
the Sun is compared with its synthetic spectrum based on 
the Holweger-Miiller model (1974) and computed by using the 
atomic line list of Sauval (2002, priv. comm.). Panel (b) shows 
the same comparison as panel (a) , but now at a resolving power 
of 1000 for the wavelength range going from 3.52 - 4.08 fj.m 
(band IE). The data of band IE of the ISO-SWS spectrum of 
a Cen A (at a resolving power of 1000) are plotted in panel 
(c). 



ery series, the fine structure is also included for the core 
calculations. The problematic computation of the self 
broadening of hydrogen lines (Barklem et al. 2000) and 
the unproper treatment of the convection by using the 
mixing-length theory (see, e.g., Asplund et al. (2000) 
where realistic ab-initio 3D radiative-hydrodynamical 
convection simulations have been used) may account 
for part of this discrepancy. 

Compared to the ISO-SWS data, the synthetic spectra 
of warm stars display a higher synthetic flux between 
the H5-9 and H5-8 hydrogen line (sec Fig. 8). From 
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Fig. 7. FTS-ATMOS intensity spectrum of the Sun (black) 
at a resolving power A/AA of 60000. The synthetic spectrum 
(grey) has been computed using the Holweger-Miiller model 
(Holweger & Miiller 1974) with = 1km s'^ The Bracket /3 
line is indicated by an arrow. Most of the other features are 
CO Av — 2 lines, some are atomic lines. 
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Fig. 8. Comparison between the observed ISO-SWS spectrum 
of a Car and its synthetic spectrum in the wavelength region 
from 3.20 to 4.08 fj.m. The Brackett-Q, two Pfund-lines (H5- 
8 and H5-9) and a number of Humphreys-lines are identified. 
The dashed line indicates the Humphreys ionisation edge. 



other SWS observations available in the ISO data- 
archive, we could deduce that this 'pseudo-continuum' 
starts arising for stars hotter than K2 (~ 4500 K). Since 
this effect is systematically seen for all warm stars in 
the sample, the origin of this discrepancy can not be 
a wrong multiplication factor. Note that band ID and 
band IE overlap each other over quite a large wave- 
length range, so that the multiplication factor of band 
IE can be accurately determined. Moreover such an ef- 
fect was not seen for the cooler K and M giants, the 
spectrum of which is dominated by OH features in this 
wavelength range, so we could reduce the problem as 
having an atomic origin. A scrutiny on the hydrogen 
lines shows that the high-excitation Humphreys-lines 
(from H6-18 on) — and Pfund-lines — are always cal- 
culated as too weak. Moreover, the Humphreys ion- 
isation edge occurs at 3.2823 ^m. Since the discrep- 
ancy does not appear above the limit (i.e. at shorter 
wavelengths) and is disappearing beyond the Brackett- 
a line, the conclusion is reached that a problem with 
the computation of the complex line profile of the 
crowded Humphreys hydrogen lines towards the series 
limit causes this discrepancy. The higher the transition, 
i.e. going from H6-18 to H6-93, the more the log gf- 
value should be increased in order to get a good match 
between observations and theoretical computations: in- 
creasing log gf with -1-0.3 dex gives us a good match for 
the H6-18 line in a Car, while we should increase log 
gf with +1.0 dex for the H6-16 line in a Car. We may 
conclude that a complex problem with the computa- 
tion of the pressure broadening will be on the origin of 
the discussed discrepancy. 



4. From 3.48 on, fringes at the end of band ID affect 
the ISO-SWS spectra of almost all stars in the sample. 

5. A clear discrepancy is visible at the beginning of band 
lA. For the warm stars, the H5-22 and H5-23 lines 
emerge in that part of the spectrum. An analogous dis- 
crepancy is also seen for the cool stars, though it is 
somewhat more difficult to recognise due to the pres- 
ence of many CO features (Fig. 9). Being present in 
the continuum of both warm and cool stars, this dis- 
crepancy is attributed to problems with the Relative 
Spectral Response Function (RSRF). A broad-band 
correction was already applied at the short-wavelength 
edge of band lA (Vandenbussche et al. 2001), but the 
problem seems not to be fully removed. At the band 
edges, the system response is always small. Since the 
data are divided by the RSRF, a small problem with 
the RSRF at these places may introduce a pronounced 
error at the band edge. 

6. Memory effects of the detectors make the calibration 
of band 2 for all the stars very difficult. These mem- 
ory effects are more severe for the cool stars, since the 
CO and SiO absorptions cause a steep increase (de- 
crease) in flux for the up (down) scan (Fig. 10). Band- 
border ratios are determined from the composite SEDs 
by Cohen or from the synthetic spectrum. The RSRFs 
for the sub-bands will therefore only be modelled well 
once there is a fool-proof method to correct SWS data 
for detector memory effects. Two other examples of this 
problematic RSRF modelling are visible in band 2B 
and band 2C. For almost all warm and cool stars, there 
is a 'dip' around 6 /im and a jump in flux level around 
9.3 ^m (Fig. 10). This latter problem is a residue of 
the correction for an instrumental absorption feature, 
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Fig. 9. Band lA of Vega (AO V), a Car (FO II) and a Tau 
(K5 III) are displayed. The ISO-SWS spectrum is plotted in 
black, the synthetic spectrum in grey. A clear discrepancy is 
visible at the beginning of this band (A < 2.407 ^m). 
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Fig. 10. Band 2 for Sirius (Al V) and 7 Dra (K5 III) are dis- 
played. The ISO-SWS spectrum is plotted in black, the syn- 
thetic spectrum in grey. Problems with the RSRF are indicated 
by an arrow. A coloured version of this plot is available in the 
Appendix as Fig. ??. 



documented by Vandenbussehe et al. (2001). The short- 
wavclcngth part of band 2A (where the CO Av ~ 1 ab- 
sorption starts) and of band 2C (with SiO Av — 1) are 
thus useless for passing a quantitative judgement upon 
the parameters influeneing the observed spectrum in 
this wavelength range, like Tcff, e(C), £(0), ... From 
the longer wavelength parts of these same bands, one 
will not be able to estimate directly fundamental stellar 
parameters, but these data enable us to see if there are 
no contradictions between this part of the ISO-SWS 
spectrum and the synthetic spectrum generated using 
the parameters adopted from literature or estimated 
from the ISO-SWS data in band 1. 

6.2. Cool stars: G2-M2 

1. The situation changes completely when going to the 
cool stars of the sample. While the spectrum of 
the warm stars is dominated by atomic-line features, 
molecules determine the spectral signature of the cool 
stars. A few of the — problematic — atomic features 
(see Sect. 6.1) can still be identified in these cool stars, 
e.g. the Mg-Si-Al-Ti-Fe spectral feature around 3.97 /im 
(Fig. 5) remains visible for the whole sample, even in 
/3Peg. 

2. One of the most prominent molecular features in band 
1 is the first-overtone band of carbon monoxide (CO, 



Av = 2) around 2.4 /im. Already from a Cen A on, CO 
lines emerge in band lA, although the atomic features 
are still more dominant for this star. It is striking that 
for all stars in the sample with spectral type later than 
G2, the strongest CO features (= band heads of ^^CO 
5-3, 6-4, 7-5, 8-6, and of ^^CO 4-2, 5-3, 6-4 and 7-5) are 
always predicted too strong compared to the ISO-SWS 
observation (see, e.g., band lA of 7 Dra in Fig. 11, 
where both SWS and theoretical data are rebinned to 
a resolving power of 1500, being the most conservative 
theoretical resolution for band lA (Leech et al. 2002)). 
The only exception is seen for the band head of ^^CO 
4-2 since the problem with the RSRF at the beginning 
of band lA (see point 5. in previous section) dominates 
in this wavelength range. 

In these oxygen-rich stars, the CO lines direct 
measure of the C abundance. From the present spectra, 
this carbon abundance can be estimated in two ways: 

— the strength of the CO spectral features is directly 
related to e(C); 

— the presence of many CO absorption lines causes a 
global depression in the spectrum up to ^ 2.9 //m 
(see Fig. 4 in Paper I), which may be used to de- 
termine £(C). 

Computing a synthetic spectrum with the carbon abun- 
dance determined from this last criterion, results how- 
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Fig. 11. Comparison between the ISO-SWS observation and 
the synthetic spectrum of 7 Dra at a resolving power of 1500 
(being the most conservative theoretical resolution for band lA 
Leech et al. 2002). A coloured version of this plot is available 
in the Appendix as Fig. ??. 



ever in the (strongest) computed CO spectral features 
being too strong compared to the ISO-SWS observa- 
tion, most visible in the strong bandheads (2 - 4%). 
It has to be noted that this mismatch occurs in band 
lA, where the standard deviation of the rebinned spec- 
trum is larger than for the other sub-bands (see Fig. 
6 in Paper I) and that the error is within the quoted 
accuracy of ISO-SWS in band lA (Leech et al. 2002). 
However, this mismatch is not random, in the sense 
that the observed CO features are always weaker than 
the synthetic ones. 

A first step towards solving this problem was the use 
of the high-resolution FTS-KP spectrum of a Boo (see 
Section 5). As was explained in Sect. 5 a high-resolution 
synthetic spectrum was generated for a Boo and re- 
binned to a resolution of 60000. The agreement be- 
tween the observed FTS-KP and synthetic spectrum 
is extremely good! One example was already shown in 
Fig. 16, another one is depicted in Fig. 17 on page 15. 
In this latter figure, the problematic subtraction of the 
atmospheric contribution causes the spurious features 
around 2.46700 /im for the summer FTS-KP spectrum 
and around 2.46735 /im for the winter FTS-KP spec- 
trum."* When scrutinising carefully the first-overtone 
CO lines in the FTS-KP spectrum, it is obvious that 
all the i^CO 2-0 fines, and almost aU the ^^CO 3-1 fines, 
are predicted as too weak (by 1 - 2 %) and not as too 

* The colour plots showing the comparison be- 
tween FTS-KP and synthetic spectrum for the wave- 
length range 2.38 - 4.08 fj,m, may be obtained at 
www.ster.kuleuven.ac.be/" leen/artikels/IS02/FTS/. Spectral 
ranges, for which the atmospheric contribution causes too 
much trouble, are omitted. 
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strong as was suggested in previous paragraph from the 
comparison between ISO-SWS and synthetic spectra! 
Also the fundamental CO fines are predicted as being 
a few percent too weak. For a better judgement of the 
KP-SWS-synthetic correspondence, the FTS-KP spec- 
trum was rebinned to the ISO resolving power. This 
was not so straightforward due to the presence of the 
spurious features originating from the problematic sub- 
traction of the atmospheric contribution. In order to 
conserve the flux, these spurious features were replaced 
by the flux value of the synthetic spectrum. This adap- 
tation is acceptable for the following reasons: 

— most of the changes were made in regions that are 
not strongly contributing to the most prominent 
features; 

— if for a particular feature the first argument does 
not apply, the change mentioned above does still not 
introduce additional errors. The FTS-KP-spectrum 
will be locally shifted towards the ISO spectrum, 
since the synthetic predictions arc (in this partic- 
ular case) systematically in between the ISO-SWS 
and FTS-KP spectrum. Thus, if statistically rele- 
vant differences are still detected, they will always 
correspond to an underestimate of the real differ- 
ences between the ISO-SWS and FTS-KP spectra. 

Fig. 18 on page 15 shows the comparison between the 
ISO-SWS, the FTS-KP and the synthetic spectrum for 
a Boo at a resolving power of 1500. It is clear that the 
differences between the strength of the CO features in 
the ISO-SWS and FTS-KP spectrum are significant. 
This is a strong indication for problems in the calibra- 
tion process. The question now arises where this cali- 
bration problem originates from. 

Firstly, it has to be noted that the flux values are un- 
reliable in the wavelength region from 2.38 to 2.40 /im 
due to problems with the RSRF of band lA; main fea- 
tures here include the CO 2-0 P18 and the CO 2-0 P21 
lines. 

Secondly, no correlation is found with the local minima 
and maxima in the RSRF of band lA. 
Defining the spectral resolution and instrumental pro- 
file for the ISO-SWS grating spectrometers, is not 
straightforward (Lorente 1998; Lutz et al. 1999). Only 
for an AOT02 observation has the instrumental profile 
now been derived quite accurately (Lutz ct al. 1999). 
The AOTOl mode introduces an additional smoothing 
which, due to the intricacies of SWS-data acquisition, 
is different from a simple boxcar smooth. The simula- 
tion of SWS-AOTOl scans gives non-gaussian profiles, 
as can be seen in Fig. 6 by Lorente (1998). Nevertheless, 
the theoretical profile of a specd-4 observation approx- 
imates a gaussian profile very closely. Therefore, since 
the instrumental profile of an AOTOl is still not ex- 
actly known, the synthetic data were convolved with 
a gaussian with FWIIM= A/resolution. This incorrect 
gaussian instrumental profile introduces an error which 
will be most visible on the strongest lines. 
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Fig. 12. Comparison between tlie AOT06 observation (revo- 
lution 538, black) and the AOTOl speed-4 observation (grey) 
of 7 Dra. Both are rebinned to a resolving power of 1500. A 
coloured version of this plot is available in the Appendix as 
Fig. ??. 
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Fig. 13. Top: Comparison between the AOTOl speed-4 obser- 
vation of 7 Dra rebinned at 1) a resolving power of 1500 (black) 
and 2) a resolving power of 1300 (grey). Bottom: Comparison 
between the ISO-SWS AOTOl observation in band lA (black) 
and its synthetic spectrum (grey) at a resolving power of 1300. 
A coloured version of this plot is available in the Appendix as 
Fig. ??. 



Only for 7 Dra an AOT06 observation, scanning this 
wavelength range, was available in the ISO data- 
archive. The comparison between the AOTOl speed-4 
and AOT06 observation, both rebinned to a resolving 
power of 1500, is an indication that the resolving power 
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of an AOTOl speed-4 observation in band lA is lower 
than 1500 (Fig. 12). The theoretical resolving power for 
an AOTOl speed-4 observation in band lA is > 1500. 
Fig. 5 in Lorente (1998) shows, however, a large devi- 
ation from this value, attributed to the fainter contin- 
uum of the source used for measuring the instrumental 
profile and the less accurate fitting in this band. The 
use of a resolving power of 1300 — instead of the theo- 
retical resolving power of 1500 — for band lA yields a) 
almost no difference from the AOTOl SWS data at a 
resolving power of 1500 and b) a better match between 
the SWS and synthetic data for the strongest CO fea- 
tures (Fig. 13). This observational resolving power of 
1300 is also in good agreement with the observational 
value given by Lorente (1998) in her Fig. 5. 
The incorrect use of a gaussian instrumental profile, 
together with too high a — theoretical — resolving 
power of 1500 form the origin of the discrepancy seen 
for the strongest CO features. Therefore, the resolving 
power of band 1 A was taken to be 1300 instead of 1500. 
Part of the other discrepancies seen in band lA may 
be explained by problems with 1. the accuracy of the 
oscillator strengths of atomic transitions in the near- 
IR (see point 1. in this section) or 2. the RSRF in the 
beginning of band lA (see point 5. in previous section). 

wavenuriber [cm ^] 
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Fig. 14. Summer and winter FTS-KP spectra of a Boo at 
a resolving power A/AA of 60000. They are compared with 
the synthetic spectrum for a Boo with stellar parameters Tcff 
= 4350 K, log g = 1.50, M = 0.75 M©, [Fe/H] = -0.50, Ct 
= 1.7km s-i, i2c/"C = 7, e(C) = 7.96, e(N) = 7.61, e(0) 
= 8.68, £(Mg) = 7.33 and Ft = 3km s~\ The OH 1-0 lines 
are predicted as too weak. A coloured version of this plot is 
available in the Appendix as Fig. ??. 



3. Concentrating in the OH-lines, we see that for both 
the high-resolution FTS-KP (Fig. 14, Fig. 16) and 
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the medium-resolution SWS spectra (Fig. 19), the 
strongest lines (OH 1-0 and OH 2-1 lines) are pre- 
dicted as too weak, while the other OH lines match 
very well. Since the same effect occurs for these two 
different observations, it is plausible to assume that 
the origin of the problem is situated in the theoreti- 
cal model or in the synthetic-spectrum computation. 
Wrong oscillator strengths for the OH lines could, e.g., 
cause that kind of problems. Being based on different 
electric dipole moment functions (EDMFs), the OH- 
line lists of Sauval (Melen et al. 1995), Partridge & 
Schwenke (1997) and Goldman et al. (1998), do how- 
ever all show the same trend. The difference in \oggf- 
value for the main branches is small (< 0.01) be- 
tween Goldman and Sauval, and is somewhat larger 
(< 0.05) between Schwenke and Sauval. Thus, no sys- 
tematic error seems to occur in the oscillator strenghts. 
Since a similar discrepancy was also noted for the low- 
excitation CO lines, some assumptions on which the 
models are based, e.g. homogeneity, hydrostatic equi- 
librium, are cleary questionable for these stars. 

4. From 3.48 /xm on, fringes at the end of band ID affect 
the ISO-SWS spectra of almost all stars in the sample. 

5. Also for the cool stars, the same remark as given in 
point 6 for the warm stars concerning the memory ef- 
fects applies. Clearly, the fundamental bands of CO 
and SiO are present for almost all cool stars, though 
no accurate quantitative scientific interpretation can 
be made due to these memory effects. 

7. Implications on calibration and modelling 

The results of this detailed comparison between observed 
ISO-SWS data and synthetic spectra have an impact both 
on the calibration of the ISO-SWS data and on the theo- 
retical description of stellar atmospheres. 

From the calibration point of view, a first conclusion is 
reached that the broad-band shape of the relative spectral 
response function is at the moment already quite accurate, 
although some improvements can be made at the begin- 
ning of band lA (sec Fig. 15) and band 2. Also, a fringe 
pattern is recognised at the end of band ID. Inaccuracies 
in the adopted instrumental profile used to convolve the 
synthetic spectra with, together with too high a resolving 
power, may cause the strongest CO Av = 2 band heads 
in the observed ISO-SWS spectrum to be weaker than the 
predicted line strength. 

Concerning the modelling part, problems with the con- 
struction of the theoretical model and the computation 
of the synthetic spectra are pointed out. The comparison 
between the high-resolution FTS-KP spectrum of a Boo 
and the corresponding synthetic spectrum revealed that 
the low-excitation first-overtone (and fundamental) CO 
lines and fundamental OH lines are predicted as being a 
few percent too weak, indicating that some assumptions, 
on which the models are based, are questionable for cool 
stars. This is not surprising: in the MARCS-code radia- 
tive equilibrium is assumed, also for the outermost lay- 
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Fig. 15. The observed ISO-SWS spectra of some stars with dif- 
ferent spectral types are divided by their synthetic spectrum 
calculated by using the atomic line list of Hirata & Horaguchi 
(1995). The ratio is rebinned to a resolving power of 250. Only 
at the beginning of band lA the same trend is visible in all 
stars, indicating that the broad-band shape of the relative 
spectral-response functions is already quite accurately known 
for the different sub-bands of band 1, with the only exception 
being at the beginning of band lA. The feature arising around 
3.85 ^m for the cool stars is an atomic feature discussed is 
Section 6.1. 



ers, implying that a temperature bifurcation, caused by 
e.g. effects of convection and convective overshoot with 
inhomogeneities in the upper photosphere, can not be al- 
lowed for. Consequently, the cores of e.g. CO lines — or 
the saturated CO lines in general — are not described 
with full success. Noting however the high level of accor- 
dance between observations and theoretical predictions for 
many molecular lines, we may conclude that the oscillator 
strengths for these molecular transitions arc now already 
accurate enough in order to use these lines to test some 
of the assumptions made in the mathematical stellar at- 
mosphere code: e.g. the temperature distribution can be 
disturbed in order to simulate a chromosphere, convection, 
a change in opacity, ... The complex computation of the 
hydrogen lines, together with the inaccurate atomic oscil- 
lator strengths in the infrared rendered the computation 
of the synthetic spectra for warm stars difficult. Therefore, 
one of us (J. S.) has derived empirical oscillator strengths 
from the high-resolution FTS-ATMOS spectrum of the 
Sun (for more details, we refer to Paper V of this series 
and to Sauval 2002). 

Although it was impossible to perform a detailed com- 
parison between observed and synthetic data in band 2, we 
may conclude that the continua of the synthetic spectra in 
this wavelength range are reliable, but that strong molecu- 
lar lines (e.g. band heads of CO Av = 1 and SiO Aw = 1) 
may be predicted a few percent too weak. Nevertheless, 
combining the synthetic spectra of both warm and cool 
stars, allowed us to test the recently developed method 
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Fig. 16. Summer (red) and winter (green) FTS-KP spectrum 
of a Boo at a resolving power of 60000. The differences be- 
tween the summer and winter FTS-KP spectrum of a Boo are 
mainly the result of noise resulting from the removal of the tel- 
luric spectrum. Due to the opposite heliocentric shifts of the 
spectra, the spectrum of Arcturus is fully recovered in spite of 
the telluric lines. The Arcturus spectra are compared with the 
synthetic spectrum (black) for a Boo with stellar parameters 
Teff = 4350 K, log g = 1.50, M = 0.75 Mq, [Fe/H] = -0.50, 
= 1.7km s-\ ^^C/^C = 7, £(C) = 7.96, e(N) = 7.61, e(0) = 
8.68, e(Mg) = 7.33, e(Si) = 7.20 and Ft = 3km s^V 
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Fig. 17. Summer (red) and winter (green) FTS-KP spectra of 
a Boo at a resolving power A/AA of 60000. They are compared 
with the synthetic spectrum (black) for a Boo with stellar 
parameters T^s = 4350 K, log g = 1.50, M = 0.75 M©, [Fe/H] 
= -0.50, Ct = 1.7km s-S i^C/i^C = 7, e(C) = 7.96, £(N) = 
7.61, £(0) = 8.68, £(Mg) = 7.33 and Ft = 3km s^^ 
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Fig. 18. Comparison between CO spectral features of the ISO- 
SWS spectrum of a Boo (black), the summer FTS-KP spec- 
trum (red), the winter FTS-KP spectrum (green) and the syn- 
thetic spectrum (blue) at a resolving power of 1500. 
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Fig. 19. Comparison between OH spectral features of the ISO- 
SWS spectrum of a Boo (black), the summer FTS-KP spec- 
trum (red), the winter FTS-KP spectrum (green) and the syn- 
thetic spectrum (blue) at a resolving power of 1500. The strong 
atomic feature around 3.87 /im is discussed in Sect. 6.1. 
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for memory effect correction (Kester et al. 2001) and to 
re-derive the relative spectral response function for bands 
1 and 2 for OLPIO (Vandcnbussche et al. 2001). The re- 
sults of these tests will be described in Paper V of this 
series (Decin et al., 2002, in preparation). In conjunction 
with photometric data, this same input data-set was used 
for the re-calibration of the absolute flux-level of the spec- 
tra observed with ISO-SWS (Shipman et al. 2001). In this 
way, both consistency and integrity were implemented. 
Moreover, the synthetic spectra of the standard sources of 
our sample are not only used to improve the flux calibra- 
tion of the observations taken during the nominal phase, 
but they are also an excellent tool to characterise insta- 
bilities of the SWS spectrometers during the post-helium 
mission. 

Although we have mainly concentrated on the dis- 
crepancies between the ISO-SWS and synthetic spectra 
— since this was the main task of this research — we 
would like to emphasise the good agreement between ob- 
served ISO-SWS data and theoretical spectra. The small 
discrepancies still remaining in band 1 are at the 1 - 2 % 
level for the giants, proving not only that the calibration 
of the (high-flux) sources has already reached a good level 
of accuracy, but also that the description of cool star at- 
mospheres and molecular line lists is quite accurate. 
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